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The objective of our study is to model the non-linear behavior of a near-critical liquid 
following a rapid change of the temperature and/or other thermodynamic parameters (pressure, 
external electric or gravitational field). The thermodynamic critical point is manifested by large, 
strongly correlated fluctuations of the order parameter (particle density in liquid-gas systems, 
concentration in binary solutions *" 3 ) in the critical range of scales. The largest critical length scale 
is the correlation radius r c . According to the scaling theory, r c increases as r c =r 0 e a when the non- 
dimensional distance €=(T-T C )/T C to the critical point decreases. Here, T c is the critical temperature, 
and the critical exponent a is =2/3. The relaxation time t(r c ) of fluctuations on the correlation radius 
length scale may be estimated near T c as x(r c )=t 0 e' <3v+z) , where the microscopic time t 0 is about 10 13 
sec, and the numerical value of the related critical exponent is approximately 3v+z=2. 

The normal gravity alters the nature of correlated long-range fluctuations when one reaches 
e-10' 5 , and correspondingly t(r c )-10‘ 3 seconds; this time is short when compared to the typical 
experimental time. In microgravity conditions, one may study a statistically homogeneous critical 
liquid which is orders of magnitude closer to T c than that in a ground-based experiment. Recently, 
the ZENO experiment 1 has set a benchmark €~10' 7 -10®, r c — 1 - 10 pm, and "c(r c )— 10 l - 10 3 s. This 
achievement manifests a qualitative change of the physical situation: that close to the critical point, 
the life-time of long-range fluctuations exceeds the typical experimental time-scale of seconds and 
minutes, and one is able to experimentally study the nonequilibrium, transient and steady, states of 
the material. On the other hand, in microgravity experiments conducted very close to the critical 
point a speeded-up equilibration may result in a long-living non-equilibrium state in some parts of 
the liquid. A theoretical model for the kinetics of those states is then needed to interpret the 
experimental data. 

Close to the critical point, a rapid, relatively small temperature change 6T«T C may perturb 
the thermodynamic equilibrium on many length scales. The critical fluctuations have a hierarchical 
structure, and the relaxation involves many length and time scales. Due to fragility of the near- 
critical equilibrium, the response of the system to such a perturbation is nonlinear one, and non- 
linearity plays the main role in the relaxation kinetics. 

Above the critical point, in the one-phase region, we consider the relaxation of the liquid 
following a sudden temperature change that simultaneously violates the equilibrium on many scales. 
Below T c , a non-equilibrium state may include a distribution of small scale phase droplets (clusters 
of one phase in the matrix of the other one); we consider the relaxation of such a droplet following 
a temperature change that has made the phase of the matrix stable. 

The transient near-critical states created by a rapid temperature change differ in important 
details from those studied in the equilibrium scaling theory and the dynamic renormalization group 
(DRG). The application of the DRG to the new situation assumes an extension of scaling 
hypothesis. We give physical arguments in favor of the proposed scenario of the hierarchical 
relaxation. The theory predicts new scaling laws for the large-time dependence of the temperature 
and for statistical characteristics of long-range fluctuations. These predictions may be tested by light 
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scattering experiments and by direct visualization techniques. An experimental test of the 
predictions will check the suggested scenario and the enhanced scaling hypothesis. 

The following physical situations are considered. 

1. A step away from the critical point where a large range of critical fluctuations disappears. 

2. A step towards the critical point where a new range of strongly interacting fluctuations appear. 

3. A phase droplet (nucleus) in a matrix of the other phase under varying conditions of the matrix 
stability. 

4. The nucleation of FCC (face centered cubic) clusters following a deep quench of a computer- 
generated "liquid". 

We start with a Langevin-type description of the system. The atomic scale fluctuations are 
assumed to play the role of a "thermal bath” for the larger scale degrees of freedom. At times t>T 0 
one characterizes this "thermal bath" by a time-dependent temperature T(t). The time x ch of the rapid 
temperature change 6T is large when compared to t 0 , but small when compared to the relaxation 
time on large scales r»r 0 : T(r)>T ch >T 0 . On the macroscopic time-scales, one considers a sharp 
temperature step at t=0. The mathematical model is as follows. To study the fluctuations of the 
order parameter (particle density in the liquid-gas system), one represents the order parameter field 
<t>(r,t) in the form 13 

<t>0v) = EV"; u) 
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the variable cut-off length A= l/k A is introduced in (1) following the Wilson- Kadanoff 
renormalization group (RG) method 2,3 . The probability W^^Dt^r) to find the equilibrium system 
in an element D(J> of the configuration space defines the effective Hamiltonian H A '- 2 : 
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In the definition (2), H <m,<:) is the microscopic Hamiltonian of the system, k b is the Boltzmann 
constant, F is the equilibrium free energy of the system, F A and F reg are the singular and the regular 
parts of the free energy respectively 2 . The sum in (2) is over all states of the system with <J)(r) a 
given function. The field <J)(r) defined by (1) represents only a small fraction of the degrees of 
freedom of the system, most of which belong to small length scales r<A. The effective Hamiltonian 
H,(<M is defined as an average over these small scale degrees of freedom. Up to a regular function 
of e, H*(4>,e) is the free energy of the system in the equilibrium state with a given configuration of 
the large scale field 4>(r). 

The fluctuations of the conserved order parameter <J>(r,t) are described by the Langevin 
equation of the form 


d<J> 

dt 




(31 


Here, T is the (conserving) kinetic coefficient. The extraneous random force f„,(r,t) models the 
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interaction of the order parameter with small-scale degrees of freedom of the system. The kinetic 
coefficient I\ and the effective Hamiltonian H x are dependent on the cut-off length A as described 
by the dynamic renormalization group (DRG). An essential condition of applicability of the method 
is the thermal equilibrium at length-scales R<A; the thermodynamic equilibrium at scales larger than 
A. is not required by the renormalization group method, the relaxation on these scales is described 
by the equations of the theory. In the course of relaxation, the fluctuations on scales R>A are non- 
equilibrium and exchange energy with smaller scales. One has to study separately a step towards 
the critical point (r c f in »r c in) which a new range of strongly interacting fluctuations appear, and 
a step away from the critical point (r c ,rin« r c.in) where a large range of critical fluctuations disappears. 
Here and below, the subscripts fm and in label the characteristics of the final and the initial (t<0) states 
of the system. 

For a large temperature change away from the critical point, r cfin «r c jn . Most important 
changes in the fluctuation picture are expected in the range of scales r c fin <r<r c in . In systems 
characterized by a conserved order parameter, the relaxation on small length-scales takes less time 
than on the large ones. At a given time t, the fluctuations on length-scales rcA^t) have approached 
the new equilibrium state while the larger scale fluctuations are still non-equilibrium. Here, A eq (t) 
is the equilibration length increasing with increasing time. At large t>x(r c fin ), A eq (t) exceeds the 
correlation radius r c fin of the final equilibrium state. We apply the DRG to find the effective 
Hamiltonian on large scales r>r c fin , t>x(r c fin ) in the form 
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Here, Xnn E X( e rin)~ e nn v<2 n) * s the susceptibility of the system in the final state (the isothermal 
compressibility for the liquid-gas system). The free field form of the resulting Hamiltonian means 
that at t>x(r c/in ) the temperature shift switches off the interaction in the range of scales r>r c fin . The 
relaxation of the long-range configuration to the final equilibrium state follows the laws of the free 
field kinetics. For the time-dependent average M k (t)=<|<j) k (t)| 2 > one obtains 
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The diffusion coefficient D fin is scale-independent at the length scales A»r cfin . We have also studied 
the shape relaxation of probability distribution function of long-range fluctuations. 

A rapid temperature increase requires more energy than a gradual heating. At large times, 
the excess energy is gradually released; in a system adiabatically insulated at t>t(r nn ), the energy 
release raises the temperature of the system. The scaling law for the large-time temperature "tail" 
has the form 

e (0-e A „ ~ SjT(r f/n )/r] C -, C. = (l*n)/2. (6) 

Here, r| = 0.033 is the critical exponent of the order parameter. 

Let us now consider a step toward the critical point: e in »e rin . In contrast to the above case, 
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the harmonics (J) k in the range r £ . in <l/k<r c fin strongly correlate in the final state, and the problem 
cannot be reduced to that of a single harmonic. The suggested scenario is sequential equilibration: 
at a given time t, the fluctuations on small scales up to a time-dependent equilibration length scale 
A^ft) w iU equilibrate while the harmonics <{> k with kcl/A^, will retain the initial magnitudes. These 
magnitudes are small compared to the final ones, and may be neglected. The statistics and the 
kinetics of the fluctuations on scales r<A eq (t) yield then the equilibrium scaling relations 2 - 3,4 : on a 
scale A<A eq (t) the relaxation time is t(A)~A 2 /D(A)~A 3+i/v , where z is the critical exponent for the 
viscosity 2,4 . The time dependence of A^ is controlled by the condition that t(r)<t for rcA^ft), and 
t(r)>t for r^>A. eq (t). The scale A^t) has then the equilibration time TfA^COi-t; this relation gives the 
scaling law for Vt>- 

V'>=rJf/T(r)F. P-~^ z (7) 


The theoretical values (see ref. 5 ) v=0.630 and z=0.063 give p=0.32. 

The sequential "cooling" of large scale fluctuations is accompanied by release of heat Q(t). 
The heat released on the scale A eq is transferred to smaller scales and through the small scale "heat 
bath" to the thermostat that maintains the temperature T(t). In a system that is adiabatically 
insulated during relaxation, the heat transferred from the large scales to the smaller ones will result 
in an increase in the temperature T(t). A rapid cooling results then in a temperature minimum. The 
scaling law for the large time "tail" in reduced temperature is 

e(0 « c,(T(r fi J/rl { . C_ = l/(3v+z). (8) 

The critical exponent £ =l/(3v+z) differs from £ + given by (6). Using the approximate theoretical 
values for the critical exponents (see ref. 5 ), one obtains for the liquid-gas and for the binary mixtures 
critical points £ + =0.5 17 and £ .=0.512. Both exponents are surprisingly close to each other and to 
the mean-field value 1/2. 

Below the critical temperature, there are two phases of the material, one of which is 
metastable except for the phase transition line. A non-equilibrium state may include droplets of one 
phase in the matrix of the other one. In the literature, this was studied for the case of a metastable 
matrix state, when the droplet is a nucleus of the new, stable, phase. We have considered a more 
general case when the embedding matrix may be metastable or stable. In the vicinity of a phase 
equilibrium line limited by the spinodals, the non-linearity of the equations for the two-phase system 
results in a finite width of the interface separating the droplet from the matrix. On both sides of the 
binodal, the droplet may be characterized by a critical radius; when the matrix is in the stable phase, 
the critical radius becomes negative. The nucleus then is always subcritical and decays with time. 
By changing the thermodynamic stability of the matrix, one is able to control the growth and decay 
of the nucleus. Outside the part of the thermodynamic plane separated by the spinodals of both 
phases, the relaxation of the droplet leads to the broadening of the boundary. Asymptotically, the 
profile of the order parameter tends to a Gaussian one, with the width R(t) growing according to the 
diffusion law R 2 ~t. 

The growth of small nuclei of the cubic central faced crystalline phase following a deep 
temperature quench to the liquid state spinodal was observed in a computer-generated "liquid". To 
discriminate clusters having the FCC local order, the statistical theory of pattern recognition was 


84 



applied. The method proposed allows one to compare nucleation kinetics in real and in simulated 
systems. 
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